UNCLASSIFIED 


Paper  Number 


INVESTIGATION  OF  PHYSICAL  AND  CHEMICAL  DELAY  PERIODS  OF 
DIFFERENT  FUELS  IN  THE  IGNITION  QUALITY  TESTER  (IQT) 


Ziliang  Zheng,  Tamer  Badawy,  Naeim  Henein 

Wayne  State  University 

Eric  Sattler 

US  Army  RDECOM-TARDEC 


ABSTRACT 

This  paper  investigates  the  physical  and  chemical  ignition 
delay  (ID)  periods  in  the  constant  volume  combustion 
chamber  of  the  Ignition  Quality  Tester  (IQT).  IQT  was  used  to 
determine  the  Derived  Cetane  Number  (DCN)  according  to 
ASTM  D6890-10a  standards.  The  fuels  tested  were  ultra  low 
sulfur  diesel  (ULSD),  jet  propellant-8  (JP-8),  two  synthetic 
fuels  of  Sasol  IPK  and  F-T  SPK  (S-8).  A  comparison  was 
made  between  the  DCN  and  cetane  number  (CN)  determined 
according  to  ASTM-D613  standards.  Tests  were  conducted 
under  steady  state  conditions  at  a  constant  pressure  of  2 1  bar, 
and  various  air  temperatures  ranging  from  778  K  to  848  K. 
The  rate  of  heat  release  (RHR)  was  calculated  from  the 
measured  pressure  trace  and  a  detailed  analysis  of  the  RHR 
trace  was  made  particularly  for  the  autoignition  process.  Tests 
were  conducted  to  determine  the  physical  and  chemical  delay 
periods  by  comparing  results  obtained  from  two  tests.  In  the 
first  test,  the  fuel  was  injected  into  air  according  to  ASTM 
standards.  In  the  second  test,  the  fuel  was  injected  into 
nitrogen.  The  point  at  which  the  two  resultant  pressure  traces 
separated  was  considered  to  be  the  end  of  the  physical  delay 
period.  The  effects  of  the  charge  temperature  on  the  total  ID  as 
defined  in  ASTM  D6890-10a  standards,  as  well  as  on  the 
physical  and  chemical  delays  were  determined.  It  was  noticed 
that  the  physical  delay  represented  a  significant  part  of  the 
total  ID  over  all  the  air  temperatures  covered  in  this 
investigation.  Arrhenius  plots  were  developed  to  determine  the 
apparent  activation  energy  for  each  fuel  using  different  IDs. 
The  first  was  based  on  the  total  ID  measured  according  to 
ASTM  standards.  The  second  was  the  chemical  delay 
determined  in  this  investigation.  The  activation  energy 
calculated  from  the  total  ID  showed  higher  values  for  lower 
CN  fuels  except  Sasol  IPK.  The  activation  energy  calculated 
from  the  chemical  delay  period  showed  consistency  in  the 
increase  of  the  activation  energy  with  the  drop  in  CN 
including  Sasol  IPK.  The  difference  between  the  two  findings 
could  be  explained  by  examining  the  sensitivity  of  the 
physical  delay  period  of  different  fuels  to  the  change  in  air 
temperature. 


1.  Introduction 

The  autoignition  of  fuel-air  mixtures  in  diesel  engines  has  a 
strong  impact  on  combustion,  performance,  fuel  economy  and 
engine-out  emissions.  The  autoignition  process  includes 
overlapping  physical  processes  that  lead  to  the  formation  of  an 
autoignitable  mixture  and  chemical  reactions  that  start  the 
combustion  process.  Many  efforts  have  been  made  to 
determine  the  effect  of  fuel  properties  on  each  of  these 
processes,  particularly  with  the  increased  interest  in  the  use  of 
renewable  fuels  to  reduce  the  dependence  on  petroleum  crude. 
Also,  alternative  fuels  derived  from  petroleum  crude  and 
synthetic  fuels  are  being  considered  for  use  in  military  diesel 
engines  [1].  The  separation  between  the  physical  and  chemical 
processes  has  been  known  to  be  difficult  to  ascertain  in  diesel 
engines  in  part  because  of  the  complexity  of  the  combustion 
chamber  geometry  which  affect  the  fuel  distribution,  and  the 
high  turbulence  in  the  charge  that  affects  the  processes  of 
spray  evaporation  and  mixing  with  the  air  [2,  3].  This  is  not 
the  case  in  constant  volume  vessels  where  the  charge  is 
quiescent  before  start  of  injection.  Furthermore,  in  engines  the 
compressed  air  properties  vary  during  the  autoignition  process 
which  makes  it  difficult  to  determine  the  effect  of  the  charge 
temperature  and  pressure  on  the  ignition  delay  (ID)  period. 
Many  correlations  between  the  ID  period  measured  in  engines 
and  charge  properties  were  developed  considering  the 
temperature  and  pressure  at  the  start  of  injection  [4,  5].  Other 
correlations  considered  the  arithmetic  mean  pressure  and 
temperature  [6],  or  the  integrated  mean  values  [2,  7,  8]  during 
the  ID  period.  Such  problems  are  reduced  in  constant  volume 
vessels  where  the  changes  in  the  charge  temperature  and 
pressure  during  the  ID  period  are  much  smaller  than  in 
engines.  Other  advantages  of  using  constant  volume  vessels  in 
the  study  of  the  autoignition  process  include  the  simplicity  of 
the  experimental  setup  and  the  smaller  amount  of  fuel  required 
for  the  experiment  as  compared  to  the  engine  setup  and 
operation.  Furthermore,  it  is  feasible  to  run  the  experiment  in 
the  constant  volume  vessel  by  injecting  the  fuel  into  nitrogen 
[4,  9,  10]  in  order  to  differentiate  between  the  physical  and 
chemical  processes,  which  would  be  very  difficult  and 
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expensive  in  engines.  In  a  nitrogen  environment,  the  chemical 
processes  that  occur  are  only  the  endothermic  reactions. 

The  major  parameters  which  affect  the  ID  of  conventional 
fuels  in  diesel  engines  are  mainly  the  air  temperature,  air 
pressure  [2,  4,  10]  and  the  fuel  molecular  structure  [3].  It  has 
been  found  that  the  ID  has  higher  sensitivity  to  changes  in  air 
temperature  than  changes  in  air  pressure  [2,  3].  Also,  the 
sensitivity  of  the  ID  to  changes  in  charge  temperature  and 
pressure  increases  with  the  drop  in  cetane  number  (CN)  of  the 
fuel  [3].  Furthermore,  the  ignition  quality  of  some  alternative 
fuels  does  not  always  correlate  well  with  CN  [11-13]  because 
of  the  wide  range  of  their  chemical  and  physical  properties 
including  composition,  volatility,  density  and  distillation  [9, 
14]. 

The  aim  of  this  investigation  is  to  gain  a  better 
understanding  of  the  effect  of  fuel  properties  on  the  physical 
and  chemical  parts  of  the  ignition  delay  period  at  different 
charge  temperatures.  The  experiments  were  conducted  in  the 
constant  volume  vessel  of  the  IQT  used  in  determining  the 
Derived  Cetane  Number  (DCN)  according  to  ASTM  D  6890- 
10a.  Additionally,  the  Arrhenius  plot  was  developed  to 
determine  the  apparent  activation  energy  of  each  fuel.  The 
investigation  covered  the  following  fuels:  (a)  ultra  low  sulfur 
diesel  (ULSD),  considered  as  the  baseline  fuel,  (b)  jet 
propellant-8  (JP-8),  a  military  petroleum  based  jet  fuel  [1,  13, 
15,  16],  (c)  Sasol  IPK,  coal  based  Fischer-Tropsch  catalytic 
process  jet  fuel  [16]  and  (d)  F-T  SPK  (S-8),  natural  gas  based 
Fischer-Tropsch  catalytic  process  jet  fuel  [15-18]. 

2.  ASTM  standards  for  cetane  number 

2.1  ASTM  D613-10a  to  determine  CN  using  the  CFR  engine 

The  cetane  number  of  a  fuel  is  determined  in  a  variable 
compression  ratio  Cooperative  Fuel  Research  (CFR)  engine 
equipped  with  a  movable  plug  in  a  cylindrical  pre-chamber 
under  a  specified  set  of  conditions  [19].  The  compression  ratio 
of  CFR  engine  is  adjustable  from  8:1  to  36:1  by  using  an 
external  hand- wheel.  The  ID  period  is  kept  constant  at  2.41 
milliseconds  which  is  equivalent  to  13  crank  angle  degrees 
(CAD)  at  an  engine  speed  of  900  rpm.  The  start  of  injection 
(SOI)  is  at  13  deg  before  top  dead  center  (TDC),  thus  start  of 
combustion  (SOC)  is  always  at  TDC,  with  the  intake  air  and 
cooling  water  temperatures  at  65°C  and  100°C  respectively  [3, 
9,  11].  The  injector  nozzle  opening  pressure  is  10.3±0.34 
MPa.  By  definition,  normal  cetane  (hexadecane  Ci6H34)  is 
assigned  a  CN  of  100,  whereas  hepta-methyl-nonane  (HMN) 
is  assigned  CN  of  15  which  defines  the  cetane  number  scale 
by  the  relationship  [3,  11,  19]: 

Cetane  Number  =  %n  cetane  +  0.15 (%HMN)  (1) 

This  test  method  has  many  disadvantages  such  as  the  need 
of  a  relatively  large  fuel  sample,  significant  time  consumption, 
difficulties  in  rating  low  ignition  quality  fuels  [3]  and 
reproducibility  and  repeatability  problems  [18,  19]. 


2.2  ASTM  D6890-10a  to  determine  DCN  using  the  IQT 

The  DCN  of  a  fuel  is  determined  according  to  ASTM 
D6890  standard  by  using  the  Ignition  Quality  Tester  (IQT) 

[20]  and  is  widely  used  to  rate  the  ignition  quality  of  fuels  for 
diesel  engine.  Also,  DCN  is  considered  as  one  of  the  critical 
properties  in  the  investigation  of  surrogates  for  aviation  fuels 

[21] .  The  tester  consists  of  a  constant  volume  cylindrical 
chamber,  electrically  heated  and  charged  with  high  purity 
compressed  air.  The  fuel  sample  is  filtered  prior  to  the  test  by 
using  a  syringe  fitted  with  a  disposable  5  pm  filter  and 
injected  into  a  pre-heated,  temperature-controlled  air.  A 
complete  sequence  comprises  1 5  preliminary  tests  followed  by 
32  tests.  The  values  for  the  ID  of  the  last  32  tests  are  averaged 
to  obtain  the  result.  This  test  method  covers  ID  periods  that 
range  from  3.1  to  6.5  ms  which  correspond  to  64  DCN  and  33 
DCN  respectively.  The  following  equation  is  used  to  calculate 
DCN  from  the  measured  ID. 

DCN  =  4.46  +  (186.6//D)  (2) 

For  fuels  that  have  ID  shorter  than  3.1  ms  or  longer  than 
6.5  ms,  DCN  can  be  calculated  with  less  precision  from  Eqn. 
3. 

DCN  =  83.99  x  (ID  -  1.512)"0-658  +  3.547  (3) 

The  IQT  test  is  simpler  than  the  CFR  engine  test  because 
it  requires  a  small  fuel  sample  of  approximately  100  mL  and 
the  test  time  is  approximately  20  min  on  a  fit-for-use 
instrument  [20,  22].  In  addition,  the  IQT  test  has  a  higher 
repeatability  and  reproducibility  as  compared  to  the  CFR 
engine  method  [23,  24],  as  shown  in  Appendix  A. 
Furthermore,  the  IQT  test  is  fully  automated  and  can  be  run 
with  minimal  operator  intervention. 

The  IQT  equipment  was  utilized  in  the  present 
investigation  as  a  platform  to  investigate  the  physical  and 
chemical  processes  that  lead  to  the  autoignition  of  different 
fuels  under  well-controlled  charge  pressures  and 
temperatures.  There  are  many  definitions  of  the  ignition  delay 
period  in  the  literature.  All  researchers  agree  on  the  start  of 
fuel  injection  as  the  start  of  ID.  However,  several  criteria  have 
been  used  to  define  the  end  of  ID  or  the  start  of  combustion 
[2-4,  7-9,  25-27].  In  the  current  investigation  the  end  of  the  ID 
is  considered  to  be  the  recovery  point  in  the  pressure  trace  as 
defined  by  ASTM  D6890-10a  standard. 

3.  IQT  setup  and  procedure 

The  IQT  is  a  bench-scale  device,  comprising  of  a  constant 
volume  combustion  chamber  with  external  electrical  heating 
elements,  a  fuel  injection  system,  an  intake  system,  an  exhaust 
system,  a  data  acquisition  system,  and  a  cooling  system.  The 
combustion  chamber  is  a  cavity  along  a  central  axis  of  the 
body,  having  a  volume  of  0.213±0.002L,  which  is  pre-heated 
to  the  standard  test  temperature  of  approximately  828  K  by 
nine  cartridge-type  resistance  heaters.  The  charge  air  pressure 
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and  temperature  are  2.137±0.007  MPa  and  818  ±  30  K 
respectively  [20].  The  IQT  is  a  high  temperature  combustion 
device  that  simulates  compression  temperatures  in  diesel 
engines,  but  at  a  reduced  air  pressure. 

A  pneumatically  driven  mechanical  fuel  pump  is  used  to 
compress  the  fuel  and  deliver  it  into  the  chamber  through  an 
inward  opening  single-hole  pintle -type  injector  nozzle.  The 
needle  is  spring-loaded  with  a  screw  and  a  lock  nut  for 
adjusting  the  nozzle  opening  pressure/release  setting  [28].  The 
air  pressure  of  the  fuel  pump  actuator  is  1.21±0.03  MPa.  The 
fuel  injection  pressure  is  approximately  22.5  MPa  during  the 
main  injection  period.  A  fixed  volume  of  fuel  is  injected.  For 
n-heptane,  the  amount  of  fuel  injected  is  72±7  mg  for  each 
injection  event  [20,  29].  A  liquid-cooled  piezo-electric 
pressure  transducer,  installed  along  the  axis  of  the  combustion 
chamber  opposite  to  the  nozzle,  measures  the  chamber  gas 
pressure  before  the  start  of  injection  till  the  end  of  the  test.  The 
definition  of  ID  for  IQT  is  the  time  elapsed  between  SOI  and 
the  SOC  defined  as  “combustion  recovery  point”  [9,  10,  14], 
as  shown  in  Figure  1,  for  n  heptane.  A  schematic  of  IQT  along 
with  a  list  of  the  other  parameters  and  running  conditions  are 
given  in  Figure  2. 

The  IQT  calibration  procedure  includes  two  reference 
fuels:  n-heptane  with  minimum  purity  of  99.5%  on  volume 
basis  and  methylcyclohexane  with  a  minimum  purity  of  99.0% 
on  volume  basis.  The  average  of  three  acceptable  ID  results  of 
heptane  is  required  to  be  within  3.78±0.01  ms,  with  no 
individual  results  below  3.72  or  above  3.84  ms.  The  average 
of  two  acceptable  ID  results  of  methylcyclohexane  is  required 
to  be  within  10.4±0.5  ms,  with  no  individual  results  below  9.8 
or  above  11.0  ms.  If  any  of  the  results  is  outside  the  specified 
range,  the  system  diagnostic  procedures  should  be  used  to 
determine  the  problem  and  a  new  calibration  is  performed  [20, 
22]. 


Figure  1  -  A  sample  of  traces  for  the  needlelift  and  chamber  pressure 
depicting  the  definition  of  IQT  ignition  delay  time 


Figure  2  -  Schematic  of  the  Ignition  Quality  Tester  setup 
at  Wayne  State  University 

4.  Results  and  Discussion 

4.1.  Comparison  between  CN  and  DCN for  different  fuels 

The  CN  and  DCN  for  different  fuels  were  measured 
according  to  ASTM  D613-10a  and  D6890-10a  standards 
respectively.  For  this  investigation,  CN  was  measured  in  the 
CFR  engine  at  Southwest  Research  Institute  (SWRI)  and  DCN 
was  measured  at  Wayne  State  University  using  the  IQT 
equipment.  The  charge  air  pressure  and  temperature  in  the 
IQT  were  kept  constant  at  21.37±0.07  bar  and  828  K 
respectively,  and  the  equipment  was  calibrated  using  n- 
heptane  fuel.  The  ID  was  measured  for  ULSD,  S-8,  JP-8  and 
Sasol  IPK  and  their  DCNs  were  subsequently  calculated.  DCN 
for  Sasol  IPK  was  calculated  by  using  Eqn.  3  as  explained  in 
the  previous  section,  because  its  ID  was  longer  than  6.5 
milliseconds.  The  DCNs  for  the  rest  of  fuels  were  calculated 
by  Eqn.  2.  The  values  for  the  CN,  DCN  and  the  ID  period  for 
the  different  fuels  are  given  in  Table  1.  Figure  3  shows  CN 
and  DCN  plotted  against  the  ID  period  in  milliseconds.  It  is 
noticed  that  CN  and  DCN  are  the  same  for  ULSD  where  the 
ID  is  5  ms.  Minor  differences  between  CN  and  DCN  are 
observed  for  S-8  and  JP-8  and  can  be  attributed  to 
inaccuracies  while  running  the  tests  on  the  CFR  engine  [18, 
19].  This  is  not  the  case  for  Sasol  IPK  where  a  relatively  larger 
discrepancy  between  the  two  methods  is  observed.  This 
discrepancy  can  be  attributed  to  the  large  uncertainty  in  Eqn.  3 
in  determining  the  DCN  for  low  CN  fuels  as  described  earlier 
[20].  Another  factor  that  contributes  to  the  large  discrepancy 
for  low  ignition  quality  fuels  is  related  to  the  differences 
between  the  engine  and  the  IQT  environment.  In  the  engine 
environment  the  gas  pressure  and  turbulence  are  much  higher 
than  in  the  stagnant  environment  of  the  IQT  [3]. 


Table  1  -  CFR  CN  and  IQT  DCN  for  different  fuels 


Fuel 

IQT  ID 

IQT  DCN 

CFR  CN 

ULSD 

4.95 

42.1 

42.3 

S-8 

3.46 

58.4 

61.0 

JP-8 

4.09 

50.1 

49.0 

Sasol  IPK 

6.94 

31.1 

25.4 
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Figure  3  -  ASTM  D613  CN  and  D6890  DCN  versus  ID  period  for 
different  fuels 


4.2.  Contributions  of  the  physical  and  chemical  processes  in 

the  autoignition  of  different  fuels 

In  order  to  differentiate  between  the  contributions  of  the 
physical  and  chemical  processes  during  the  autoignition  of 
different  fuels,  tests  were  conducted  in  two  steps.  The  first 
was  with  fuel  injection  into  air  and  the  second  was  with  fuel 
injection  into  nitrogen,  under  the  conditions  specified  by 
ASTM  D6890  standard.  For  the  nitrogen  case,  the  physical 
processes  and  some  endothermic  reactions  contribute  to  the 
rate  of  heat  release.  The  differences  in  the  pressure  between 
the  two  tests  provide  information  regarding  the  start  of  the 
active  exothermic  reactions  leading  to  autoignition  and 
combustion.  In  both  cases,  the  rate  of  heat  release  (RHR)  was 
computed  from  the  pressure  trace  [30].  Figure  4  shows  the 
details  of  the  needlelift  (N.L),  pressure  trace  and  RHR  during 
the  ID  period.  The  oscillations  in  the  pressure  and  RHR  traces 
make  it  difficult  to  determine  the  start  of  the  exothermic 
reactions. 


A  low  pass  filter  was  applied  to  the  raw  pressure  signal 
with  a  cutoff  frequency  of  6  kHz  to  eliminate  the  noise  as 
explained  in  Appendix  B.  Further,  a  band  pass  filter  was 
applied  to  the  pressure  and  N.L  traces  for  the  four  fuels  with 
low  and  high  cutoff  frequencies  of  one  and  four  kHz 
respectively  as  shown  in  Figure  5a.  The  figure  shows 
similarity  between  the  oscillations  in  the  pressure  and  the  N.L 
traces,  but  with  a  slight  phase  shift.  It  is  noticed  that  the 
oscillations  are  the  same  for  all  the  fuels  in  spite  of  the 
differences  in  their  physical  properties,  given  in  Table  A.l  in 
Appendix  A.  Figure  5b  represents  a  spectral  analysis  of  the 
filtered  pressure  and  N.L  signals  [31].  It  is  clear  that  the  two 
traces  have  similar  spectra  that  peak  at  a  frequency  of  2.2  kHz. 
This  indicates  that  the  oscillations  in  the  pressure  trace  are 
related  to  the  injection  process.  This  relationship  was 
investigated  in  National  Renewable  Energy  Laboratory 
(NREL)  [29,  32]  by  monitoring  the  N.L  and  the  corresponding 
changes  in  high  speed  images  taken  for  the  fuel  spray  during 
the  injection  process.  NREL  concluded  that  the  oscillations 
were  caused  by  the  changes  in  the  rate  of  fuel  delivery  during 
the  injection  process.  Also,  the  changes  in  the  rate  of  fuel 
delivery  were  caused  by  changes  in  the  flow  area  in  the  pintle 
nozzle,  in  addition  to  oscillations  observed  in  the  needle 
position  during  the  injection  process.  It  is  worth  noting  that  a 
detailed  analysis  of  the  pressure  traces  for  the  32  tests,  used  to 
determine  the  DCN  in  this  investigation,  indicated  test-to-test 
variability  in  these  oscillations.  Such  variability  can  be 
explained  by  the  test-to-test  variability  in  the  needle 
oscillations  reported  by  NREL  [29,  32]. 

In  the  present  investigation,  the  oscillations  in  the  pressure 
trace  were  filtered  in  order  to  differentiate  between  physical 
and  chemical  ignition  delay  periods  as  explained  in  Appendix 
B. 


Figure  5a  -Filtered  pressure  and  needlelift  signals  for  four  fuels 


Figure  4  -  Raw  data  of  N.  L.,  pressure  and  RHR  traces  for  n-heptane 
4.2.1.  Analysis  of  the  oscillations  in  the  pressure  trace 
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Figure  5b  -  Spectral  analysis  for  the  filtered  signals  of  pressure  and 
needlelift  for  different  fuels 

4.2.2.  Measurement  of  Physical  and  Chemical  Ignition 
delays 

The  physical  and  chemical  ignition  delays  were 
determined  by  comparing  the  filtered  pressure  and  RHR  traces 
for  fuel  injection  in  air  with  the  corresponding  traces  for  fuel 
injection  in  nitrogen.  Figure  6  shows  filtered  traces  for 
pressure,  RHR,  gas  average  temperature  and  N.L  for  injection 
into  air  and  nitrogen.  The  N.L  traces  for  both  cases  show  the 
SOI  at  1.8  ms.  The  drop  in  pressure  immediately  following  the 
SOI  is  a  measure  of  the  drop  in  temperature  caused  by  the  heat 
transfer  from  the  charge  to  the  fuel  for  evaporation  and 
endothermic  reactions.  It  is  observed  that  both  the  gas  pressure 
traces  and  RHR  traces  for  the  injection  in  air  and  in  nitrogen 
coincide  during  the  period  from  SOI  to  4.25  ms,  after  which 
they  separate.  The  trace  for  injection  into  nitrogen  continues  to 
drop  indicating  a  continuation  of  liquid  evaporation  and 
endothermic  reactions.  The  trace  for  injection  into  air  changes 
its  slope  indicating  the  start  of  active  exothermic  reactions. 
The  point  of  separation  between  the  two  traces  is  defined  as 
point  of  inflection  (POI)  [7]  determining  the  start  of  active 
exothermic  reactions  and  the  end  of  the  physical  delay  period 
[9,  10].  While  the  exothermic  reactions  start  after  POI,  the 
trace  for  injection  into  air  continues  to  drop  and  reaches  a 
valley  denoted  by  ‘b\  This  indicates  that  endothermic 
reactions  are  predominant  between  POI  and  point  ‘b\  After 
point  4b’,  the  trace  starts  to  rise  up,  indicating  the  domination 
of  the  exothermic  reactions  that  lead  to  combustion.  It  should 
be  noted  that  the  contribution  of  the  physical  delay  in  the  total 
ignition  delay  is  more  significant  in  the  IQT  compared  to 
diesel  engines  [33].  This  is  mainly  caused  by  the  differences  in 
the  injection  process,  the  charge  pressure  and  turbulence  in 
two  experimental  set  setups. 
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Figure  6  -  Traces  for  needlelift,  pressure,  RHR  and  temperature  for 
injection  of  heptane  in  air  and  in  nitrogen 

4.2.3.  Effect  of  fuel  properties  on  the  physical  and  chemical 

ignition  delays  in  the  IQT  environment 
The  time  histories  of  the  gas  pressure,  N.L,  gas  mass 
average  temperature  and  RHR  are  given  in  Figure  7  for  four 
fuels.  In  this  analysis,  ULSD  is  considered  to  be  the  baseline 
fuel  for  the  comparison  with  other  fuels.  The  figure  shows 
that,  for  all  the  fuels,  SOI  occurs  at  1.8  ms  and  the  main 
injection  process  is  completed  in  about  2  ms  before  the  rise  in 
pressure  due  to  combustion.  As  expected,  S-8  has  the  shortest 
ID  because  it  has  the  highest  volatility  and  DCN.  The 
volatility  of  a  fuel  is  inversely  proportional  to  its  flash  point 
[34-36].  Table  A.l  gives  the  flash  point  for  the  fuels.  The  flash 
point  of  S-8  is  44.5°C  compared  to  75°C  for  ULSD.  High 
volatility  enhances  fuel  evaporation  and  the  formation  of  an 
ignitable  mixture  [7,  37],  and  shortens  the  physical  delay 
period.  Again,  the  volatility  and  DCN  contribute  to  the  shorter 
ID  of  JP-8  compared  to  ULSD.  The  flash  point  of  JP-8  is 
49.5°C  compared  to  75°C  for  ULSD.  The  DCN  for  JP-8  is 
50.1  compared  to  42.1  for  ULSD.  Sasol  IPK  has  the  lowest 
DCN  because  of  its  longest  ID  period  compared  to  the  other 
fuels,  in  spite  of  its  higher  volatility,  compared  to  ULSD  and 
JP-8.  The  drop  in  the  mass  average  temperature  caused  by 
energy  consumed  by  fuel  evaporation  and  endothermic 
reactions  is  more  evident  for  Sasol  IPK  compared  to  the  other 
fuels.  Figure  7  shows  that  the  rise  in  pressure  due  to  the 
combustion,  which  is  mostly  in  the  premixed  mode,  is  in  line 
with  the  length  of  the  ID  for  the  different  fuels.  Figure  8  gives 
the  details  of  RHR  during  the  ID  period.  It  is  observed  that  S- 
8,  ULSD  and  JP-8  depict  the  one-stage  autoignition  process, 
while  Sasol  IPK  depicts  the  two-stage  autoignition.  The 
pressure  rise  due  to  combustion  depends  mainly  on  the  length 
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of  the  ID  and  volatility  the  fuel.  S-8  has  the  shortest  ID  and 
produced  the  lowest  peak  pressure  among  all  fuels.  Also,  the 
pressure  keeps  increasing  at  a  slow  rate  indicating  diffusion 
combustion.  Figure  9  shows  the  normalized  RHR  integral  for 
different  fuels,  from  the  start  of  injection  to  the  end  of 
combustion.  Table  2  gives  the  heat  released  in  premixed  and 
diffusion  combustion  and  the  premixed  combustion  fraction. 
S-8  has  the  lowest  premixed  combustion  fraction  among  all 
the  fuels,  while  Sasol  IPK  has  the  highest  fraction.  The  sharp 
drop  in  the  integrated  RHR  for  Sasol  IPK  is  the  result  of  the 
release  of  almost  all  the  energy  in  the  premixed  combustion 
mode. 


Time  (ms) 


Figure  7  -  Pressure,  needlelift,  RHR,  temperature  traces  for  four 
fuels  tested  at  temperature  828  K 


Figure  8  -  Details  of  RHR  traces  during  autoignition  for  four  fuels  at 
temperature  828  K 


Figure  9  -  Normalized  integral  RHR  of  four  fuels 
at  temperature  82 8K 


Table  2  -  Energy  Release  of  four  fuels  at  temperature  828  K 


Fuel 

Pre-mixed  (kJ) 

Diffusion  (kJ) 

Pre-mixed  fraction 

ULSD 

2.2506 

0.6582 

77.37% 

S-8 

1.6499 

0.7757 

68.02% 

JP-8 

1.9135 

0.8107 

70.24% 

Sasol  IPK 

2.7904 

0.1393 

95.25% 

4.2.4.  Effect  of  charge  temperature  on  the  physical  and 
chemical  ignition  delays  of  different  fuels 
A  series  of  tests  were  conducted,  at  charge  temperatures 
between  778K  and  848  K  by  controlling  the  chamber  skin 
temperature.  All  the  tests  were  conducted  at  a  constant  charge 
pressure  of  21.37±0.07  bar  (310±1  psi).  The  traces  were 
analyzed  and  the  physical  and  chemical  delays  were 
determined  for  each  of  the  fuels  at  different  temperatures  and 
the  results  are  given  in  Table  3.  The  data  are  plotted  in  Figure 
10  versus  the  integrated  mean  temperature  during  the  ignition 
delay  period,  rather  than  the  temperature  at  the  start  of 
injection  [2]. 
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Figure.lOa  -  Physical  ID  vs  integrated  mean  temperature  during  the 
physical  process 

Figure  10b  -  Chemical  ID  vs  integrated  mean  temperature  during  the 
chemical  process 


Table  3  -  Physical  and  chemical  delay  for  four  fuels  at  different  temperature 


ULSD 

S-8 

JP-8 

Sasol  IPK 

Test  T  (K) 

Phy  ID 

(ms) 

Chem  ID 
(ms) 

Phy  ID 

(ms) 

Chem  ID 
(ms) 

Phy  ID 

(ms) 

Chem  ID 
(ms) 

Phy  ID 

(ms) 

Chem  ID 
(ms) 

778 

4.8 

2.12 

3.38 

1.48 

3.84 

1.9 

5.58 

3.1 

788 

4.26 

2.1 

3.06 

1.32 

3.58 

1.72 

4.92 

2.8 

798 

4 

1.86 

2.94 

1.26 

3.38 

1.52 

4.76 

2.68 

808 

3.68 

1.8 

2.58 

1.2 

3.16 

1.52 

4.5 

2.5 

818 

3.58 

1.64 

2.5 

1.1 

2.94 

1.38 

4.22 

2.44 

828 

3.34 

1.56 

2.36 

1.1 

2.66 

1.38 

4.06 

2.34 

838 

3.08 

1.52 

2.22 

1.1 

2.56 

1.34 

4.02 

2.06 

848 

2.98 

1.42 

2.08 

1.06 

2.46 

1.24 

3.96 

2 

Reduced 

Percentage 

38% 

33% 

38% 

28% 

36% 

35% 

29% 

36% 

Figure  1 1  shows  the  chamber  pressure  and  temperatures 
before  the  start  of  injection  and  the  integrated  mean 
temperature  during  the  ID  delay  period.  It  should  be  noted  that 
the  air  temperature  along  the  fuel  spray  should  be  lower  than 
the  mass  average  temperature,  calculated  from  the  pressure 
trace  and  the  mass  of  the  charge. 


Figure  12  is  a  3D  plot  of  the  total,  physical  and  chemical  ID 
periods  plotted  versus  the  integrated  mean  temperature  for 
different  DCN.  It  is  observed  that  the  physical  delay  is  a  major 
part  of  the  total  ID  period  for  all  the  fuels  in  the  IQT.  It 
indicates  that  the  DCN  does  not  only  represent  the  chemical 


T  mean  (K) 


Figure.ll  -  Chamber  pressure  &  temperatures  before  SOI  and  mean 
pressure  vs  integrated  mean  temperature  during  the  ID  period 


characteristics  of  the  fuel  but  also  takes  into  consideration  the 


physical  properties  of  the  fuel.  Further,  it  is  observed  that  the 
physical  delay  is  more  sensitive  to  the  increase  in  temperature, 
particularly  for  the  higher  DCN  fuels,  as  discussed  in  detail  in 
the  next  section. 
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Figure.12  -  3D  plot  of  the  ID  period,  integrated  mean  temperature 
and  DCN 

4.3.  Forms  of  Arrhenius  plot  and  corresponding  apparent 
activation  energy  for  different  fuels 

In  this  section  different  forms  of  the  Arrhenius  plot  are 
demonstrated  and  the  corresponding  apparent  activation 
energies,  Ea,  are  calculated  for  different  fuels.  Ea  is 
determined  from  Eqn.  4  and  its  value  depends  on  the 
temperature,  T,  used  in  the  equation.  Temperatures  of 
significance  can  be  at  the  SOI  [4,  5],  the  arithmetic  mean 
temperature  during  ID  [6]  and  the  integrated  mean 
temperature  during  the  ID  [2,  7,  8].  The  first  form  of 
Arrhenius  plot  is  for  the  total  ID  versus  the  reciprocal  of  the 
charge  temperature  before  the  start  of  injection  and  integrated 
mean  temperature  during  the  total  ID  period,  as  Figure  13  a  & 
b..  The  results  of  this  form  would  be  of  use  for  engine  design 
where  diesel  cycle  simulation  is  used  for  determining  overall 
engine  performance.  Also,  this  information  can  be  used  for 
electronic  engine  control,  where  fast  computations  are  needed 
for  ECU  operation.  The  second  plot  is  for  the  chemical  ID 
versus  the  reciprocal  of  the  integrated  mean  temperature 
during  its  ID  period,  as  Figure  14.  Here  the  results  would  be  of 
value  in  the  validation  of  detailed  combustion  mechanisms 
developed  for  different  fuels.  Also,  such  information  can  be 
used  for  the  development  of  surrogates  for  different  fuels. 
Since  the  chemical  reactions  are  very  sensitive  to 
temperatures,  it  is  felt  that  the  temperature  between  the  start  of 
injection  and  the  end  of  ID  would  be  the  most  appropriate. 
The  third  plot  is  for  the  physical  delay  versus  the  integrated 
mean  temperature  during  its  ID  period,  as  Figure  15.  This  plot 
explains  the  discrepancy  between  the  apparent  activation 
energy  calculated  from  the  first  and  the  second  plots. 

Ln(ID )  =  LnA  H -  (4) 

Ru  T 

4.3.1.  Apparent  activation  energy  for  the  total  ID  plotted 
versus  the  charge  temperature  before  the  start  of 


injection  and  integrated  mean  temperature  during  the 
total  ID  period 

Figure  13  a  &  b  are  the  plot  of  the  total  ID  measured  from 
the  start  of  injection  to  the  pressure  recovery  point,  specified 
in  ASTM  Standard  6890- 10a,  as  end  of  ID.  The  apparent 
activation  energies  calculated  from  Figure  13  are  given  in 
Table  4.  It  is  noticed  that  for  all  the  fuels  except  Sasol  IPK,  Ea 
increases  with  the  drop  in  DCN.  This  is  the  trend  in  many 
chemical  reactions.  However,  this  trend  is  not  followed  by 
Sasol  IPK.  Similar  observations  have  been  reported  in 
investigations  conducted  on  the  same  fuels  in  a  high  speed 
direct  injection  diesel  engine  [7]  and  on  the  FACE  (fuels  for 
advanced  combustion  engines)  diesels  in  IQT  by  NREL  [29]. 


Table  4  -  Apparent  activation  energy  based  on  total  ignition  delay 
for  different  fuels 


Fuel 

Ea  (kJ/mole)  from  Tcharge 

Ea  (kJ/mole)  from  Tmean 

ULSD 

15.81 

28.56 

S-8 

14.75 

27.56 

JP-8 

15.52 

28.66 

Sasol  IPK 

11.65 

22.14 

Figure  13a  -  Arrhenius  plots  for  total  ignition  delay  versus  the 
charge  temperature  before  SOI  for  different  fuels 
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Figure  13b  -  Arrhenius  plots  for  total  ignition  delay  versus  the 
integrated  mean  temperature  during  ID  period  for  different  fuels 

4.3.2.  Apparent  Activation  energy  for  the  chemical  ID 
plotted  versus  the  integrated  mean  charge 
temperature  during  the  chemical  ID 
Figure  14  shows  a  plot  of  the  chemical  ID  versus  the 
integrated  mean  charge  temperature  during  the  chemical  ID.  It 
should  be  noted  that  the  apparent  activation  energies 
calculated  from  Figure  14  are  not  for  elementary  chemical 
reactions  but  rather  for  the  global  of  a  large  number  of  known, 
unknown,  simple  and  complex  chemical  reactions  that  occur 
during  the  autoignition  process  within  the  temperature  range 
studied  [2].  The  apparent  activation  energies,  Ea,  of  different 
fuels  are  calculated,  as  shown  in  Table  5. 


Table  5  -  Apparent  activation  energy  based  on  chemical  ignition 
delay  for  different  fuels 


Fuel 

Apparent  activation  energy 

Ea  (kJ/mole) 

ULSD 

27.31 

S-8 

19.26 

JP-8 

24.77 

Sasol  IPK 

28.57 

4.3.3.  Apparent  Activation  energy  for  the  physical  ID 
plotted  versus  the  integrated  mean  charge 
temperature  during  the  physical  ID 
Figure  15  shows  the  data  in  an  Arrhenius  plot  for  the 
physical  ignition  delay  given  in  Figure  10.  The  slopes  of  the 
traces  in  Figure  15  are  given  in  Table  6  and  show  that  the 
physical  processes  of  higher  DCN  fuels  are  more  sensitive  to 
changes  in  the  charge  temperature  than  the  fuels  with  a  lower 
DCN.  This  can  be  explained  by  the  fact  that  the  start  of  the 
exothermic  reactions  signifies  the  end  of  physical  ID.  In 
general,  chemical  reactions  are  more  sensitive  to  increase  in 
temperature  which  makes  the  physical  ID  equally  sensitive  to 
the  increase  in  temperature. 
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Figure  15  -  Arrhenius  plots  for  physical  ignition  delay  versus  the 
integrated  mean  temperature  for  different  fuels 

Table  6  -  Apparent  activation  energy  based  on  physical  ignition 
delay  for  different  fuels 


Apparent  activation  energy 
Ea  (kJ/mole) 


ULSD 

30.60 

S-8 

33.43 

JP-8 

32.01 

Sasol  IPK 

23.31 

Figure  14  -  Arrhenius  plots  for  chemical  ignition  delay  versus  the 
integrated  mean  temperature  for  different  fuels 

Table  5  shows  that  S-8  which  has  the  shortest  chemical 
ID  and  the  highest  DCN  has  the  lowest  activation  energy, 
19.26  kJ/mole.  Meanwhile,  Sasol  IPK  which  has  the  longest 
ID  and  lowest  DCN  has  the  highest  activation  energy,  28.57 
kJ/mole.  The  activation  energy  of  ULSD  is  27.31  kJ/mole 
which  is  higher  than  24.77  kJ/mole  for  JP-8. 

The  question  here  is  why  there  is  a  difference  between  the 
trends  shown  in  Figure  13  and  14,  particularly  for  low  DCN 
fuel.  One  of  the  factors  that  might  have  contributed  to  this 
difference  is  the  sensitivity  of  the  physical  processes  to  the 
change  in  charge  temperature. 


4.3.4.  Arrhenius  parameters  of  the  total,  physical  and 
chemical  ID  for  different  fuels 
By  using  the  ID  and  temperature  during  the  different 
process  period,  the  Arrhenius  parameters  of  total,  physical  and 
chemical  ID  period  are  determined  respectively  for  different 
fuels  using  Eqn.  4,  as  shown  in  Table  7.  It  shows  that  the  pre¬ 
exponential  factors  (A)  and  activation  energies  are  different 
for  the  total  ID,  as  the  temperature  at  SOI  and  during  the  total 
ID  period  are  taken  into  the  calculations  respectively.  It 
indicates  that  the  temperature  should  be  chosen  appropriately 
in  order  to  determine  the  apparent  activation  energy  of  a 
certain  interest. 

Since  the  IQT  is  a  heterogeneous  combustion  system,  the 
measured  total  ID  includes  both  physical  process  and  chemical 
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reactions.  As  a  result,  the  total  ID  can  be  expressed  in  Arrhenius  equation  as  Eqn.5. 

EVhy  Echem 

ID  total  ~  IDphy  +  IDchem  —  Aphy  '  e  u  mean-Phy  +  Acfiem  ■  e  u  mean_chem  (5) 

Table  7  -  Arrhenius  parameters  of  ignition  delay  for  different  fuels 


ULSD  S-8  JP-8  Sasol  IPK 

Temperature  (K)  A  E  (kJ/mol)  A  E  (kJ/mol)  A  E  (kJ/mol)  A  E  (kJ/mol) 

Tsoi  0.1628  15.81  0.1438  14.75  0.1459  15.52  0.5660  11.65 

Tmean  total  0.0674  28.56  0.0577  27.56  0.0578  28.66  0.2432  22.14 

Tmeanphy  0.0346  30.60  0.0167  33.43  0.0238  32.01  0.1232  23.31 

Tmean  chem  0-0247  27.31  0.0614  19.26  0.0321  24.77  0.0281  28.57 


5.  Conclusions 

The  following  conclusions  are  based  on  an  experimental 
investigation  of  the  role  of  the  physical  and  chemical 
processes  in  the  autoignition  of  fuels  with  different  physical 
and  chemical  properties.  Four  fuels  were  tested  at  different 
temperatures,  varying  from  778  to  848  K.  The  physical  and 
chemical  ignition  delays  were  determined  by  two  tests.  In  the 
first  test,  the  fuel  was  injected  into  air  charge  following  the 
ASTM  D6890-10a  procedure.  In  the  second  test,  the  fuel  was 
injected  into  pure  nitrogen  under  the  same  conditions  as  the 
first  test.  The  difference  between  the  pressure  traces  measured 
in  the  two  tests  indicated  the  time  when  the  chemical 
exothermic  reactions  started  which  was  considered  to  be  the 
end  of  the  physical  delay  period.  The  total  ID  defined  by 
ASTM  D6890-10a,  the  chemical  and  physical  IDs  were 
correlated  with  the  integrated  mean  temperature  using  the 
Arrhenius  plot,  from  which  the  apparent  activation  energy  was 
determined  for  each  fuel. 

1.  A  frequency  analysis  of  the  oscillations  in  the  gas 
pressure  trace  indicated  that  they  are  related  to  the 
injection  process,  as  has  been  reported  earlier  in  the 
literature  by  using  imaging  techniques.  These  oscillations 
were  filtered  before  the  analysis  was  made. 

2.  The  low  DCN  fuel,  Sasol  IPK,  exhibited  the  two-stage 
autoignition  process,  particularly  at  the  low  charge 
temperatures,  which  was  not  the  case  for  the  other  fuels. 

3.  The  physical  delay  period  in  IQT  is  a  major  part  of  the 
total  ignition  delay  period  for  all  the  fuels.  This  is  not  the 
case  in  high  speed  direct  injection  diesel  engines  where 
the  physical  delay  period  is  a  small  fraction  of  total  delay 
period,  as  a  result  of  major  differences  in  the  design  of  the 
injection  system  and  combustion  chamber  in  the  two 
cases. 

4.  The  apparent  activation  energy  calculated  from  the  total 
ignition  delay  has  an  inverse  relationship  with 
the  DCN  of  the  fuels  for  DCNs  higher  than  40.  This  is  not 
the  case  with  Sasol  IPK  which  has  the  lowest  DCN  and 
the  lowest  apparent  activation  compared  to  other  fuels. 

5.  The  apparent  activation  energy  calculated  from  the 
chemical  delay  has  an  inverse  relationship  with  the  DCN 
of  all  fuels  including  Sasol  IPK. 


6.  The  difference  between  the  two  findings  is  explained  by 
plotting  the  physical  delay  in  an  Arrhenius  format.  The 
plots  indicated  a  higher  sensitivity  of  the  physical  delay  of 
high  DCN  fuels  to  changes  in  temperature  compared  to 
the  sensitivity  of  lower  DCN  fuels. 
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Definition/ Abbreviation 


A 

Pre-exponential  factor 

AET 

Advanced  Engine  Technology 

ARHR 

Apparent  Rate  of  Heat  Release 

CAD 

Crank  Angle  Degree 

CFR 

Cooperative  Fuel  Research 

Cl 

Compression  Ignition 

CN 

Cetane  Number 

DCN 

Derived  Cetane  Number 

Ea 

Apparent  Activation  Energy 

EGR 

Exhaust  Gas  Recirculation 

FACE 

Fuels  for  Advanced  Combustion 
Engines 

F-T  SPK 

Fischer-Tropsch  Synthetic 
Paraffinic  Kerosene 

HMN 

Hepta-methyl-nonane 

ID 

Ignition  Delay 

IQT 

Ignition  Quality  Tester 

JP-8 

Jet  Propellant  8 

LT 

Low  Temperature 

NTC 

Negative  Temperature  Coefficient 

NREL 

National  Renewable  Energy 
Laboratory 

POI 

Point  of  Inflection 

RHR 

Same  as  ARHR 

R« 

Universal  Gas  Constant 

SOI 

Start  of  Injection 

soc 

Start  of  Combustion 

SwRI 

Southwest  Research  Institute 

TDC 

Top  Dead  Center 

ULSD 

Ultra-Low- Sulfur-Diesel 
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Table  A.l  -  Fuel  Properties 


Fuel 

ULSD 

JP-8 

Sasol  IPK 

S-8 

D613  CN 

42.3 

49 

25.4 

61 

D6890  DCN 

42.1 

50.1 

31.1 

58.4 

Flash  point  (°C)  min 

75 

49.5 

42 

44.5 

Density  (@15°C)  kg/m3 

842 

798 

755 

749 

Heating  Value  (MJ/kg) 

41.2 

43.2 

44.0 

44.1 

Figure  A.l  -  Distillation  curves  for  four  fuels 


Table  A.2 

Repeatability  and  Reproducibility  Values  for  Information 


ID  (ms) 

Repeatability  (r) 

Reproducibility  (r) 

3.1 

0.030 

0.158 

3.6 

0.055 

0.198 

4.2 

0.085 

0.246 

4.8 

0.115 

0.293 

6.0 

0.175 

0.388 

6.5 

0.200 

0.428 

DCN 

Repeatability  (r) 

Reproducibility  (r) 

33 

0.67 

1.96 

40 

0.77 

2.23 

45 

0.83 

2.43 

50 

0.90 

2.62 

55 

0.96 

2.81 

60 

1.03 

3.00 

64 

1.08 

3.16 
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Appendix  B 

The  time  history  of  the  needlelift  (N.L.),  raw  pressure  and  rate  of  heat  release  (RHR)  are  shown  in  Figure  4.  It  is  observed  that  the 
pressure  trace  has  oscillations  and  signal  noise.  Moreover,  the  RHR  trace  shows  more  noise,  since  it  is  derived  from  raw  pressure 
trace.  The  technology  of  filtering  is  applied  to  eliminate  the  noise  and  oscillations,  which  enhances  the  understanding  for  the 
autoignition  period.  Figure  B.la  shows  the  zoomed-in  raw,  smoothed,  noise  filtered  and  N.L  oscillation  filtered  pressure  traces  of  n- 
heptane.  Figure  B.lb  shows  the  Fast  Fourier  Transform  (FFT)  analysis  plot  for  the  raw,  smoothed,  noise  filtered  and  N.L  oscillation 
filtered  pressure  signals.  The  solid  blue  line  is  the  raw  pressure  signal.  The  dash  red  line  is  the  smoothed  pressure  signal  by  eleven- 
point  moving  average  method.  The  dash-dot  green  line  is  the  noise  filtered  pressure  signal.  The  dash-two-dot  orange  line  is  the  N.L 
oscillation  filtered  pressure  signal.  The  frequencies  of  N.L,  combustion,  resonance  and  etc.  can  be  differentiated  and  determined  using 
wave  equation  in  the  combustion  chamber  [31,  38-42]. 

1 .  Signal  Smoothing 

The  eleven-point  moving  average  smoothing  method  is  applied  to  smooth  the  raw  pressure  as  shown  in  Figure  B.la.  From 
which,  the  SOC  can  be  determined  from  the  pressure  recovery  point.  It  is  observed  that  the  signal  noise  is  minimized  without 
any  shift  in  the  pressure  trace.  However,  all  frequencies  are  still  present  in  the  signal  even  after  smoothing  as  shown  in  Figure 
B.lb.  Further,  the  oscillations  caused  by  N.L.  are  still  observed  in  the  smoothing  trace  which  affects  the  precision  of 
determining  point  of  inflection  (POI).  In  addition,  a  very  important  disadvantage  of  the  moving  average  method  is  that  the 
point  of  the  moving  average  would  need  to  be  tuned  for  each  cycle  in  order  to  obtain  optimum  smoothing  results. 

2.  Finite  impulse  response  (FIR)  filter  with  phase  correction 

On  the  other  hand,  using  FIR  filter  with  phase  correction  is  recommended. 

A  low  pass  filter  with  a  cutoff  frequency  of  6  kHz  eliminates  the  noise  present  in  the  raw  pressure  signal  as  shown  in  Figure 
B.la.  It  is  also  observed  the  noises  with  higher  than  6  kHz  frequencies  are  wiped  out  in  the  corresponding  FFT  plot  as  shown 
in  Figure  B.lb. 

Since  the  N.L.  frequency  is  determined  at  2.2  kHz,  an  additional  low  pass  filter  with  a  cutoff  frequency  of  1.6  kHz  is  applied 
to  get  rid  of  the  N.L  oscillations  presented  in  the  pressure  trace  as  shown  in  Figure  B.la.  As  a  result,  it  is  also  observed  the 
frequencies  with  higher  than  1.6  kHz  are  eliminated  in  the  corresponding  FFT  plot  as  shown  in  Figure  Bl.b., However,  a  shift 
occurred  in  the  combustion  phase  compared  with  raw  pressure  signal  as  shown  in  Figure  B.la  and  B.lb.  Although  the  SOC 
cannot  be  accurately  determined  from  the  1.6  kHz  low  pass  filtered  signal,  the  POI  in  the  autoignition  period  can  be  defined 
with  high  accuracy. 

Thus,  the  SOC  should  be  determined  from  the  pressure  signal  after  applying  the  noise  filter  with  the  cutoff  frequency  of  6  kHz.  And 
the  POI  could  be  defined  from  the  pressure  signal  after  applying  the  N.L  oscillation  filter  with  the  cutoff  frequency  of  1.6  kHz. 
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Figure  B.la  -  Zoomed  raw,  smoothed,  noise  filtered  and  N.L  oscillation  filtered  pressure  traces  for  n-heptane 
Figure  B.lb  -  FFTs  analysis  for  the  raw,  smoothed,  noise  filtered  and  N.L  oscillation  filtered  pressure  signals 
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